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In a paper read before the Society on June 11, Sir William Thomson expressed a 
doubt as to the general truth of the Maxwell-Boltzmann doctrine concerning the 
distribution of energy among a great number of mutually acting bodies, and suggested 
that certain test cases should be investigated. The test that he proposed on that 
occasion was a number of hollow elastic spheres, each of mass M, and each containing 
a smaller elastic sphere of mass m, free to move within a larger one. This pair he 
calls a doublet. This case is within the general proof of the doctrine given below. 
It is, however, I think amenable to a simpler treatment, which has been applied to 
the case of elastic spheres external to one another. 

1. Every doublet has a centre of inertia of the sphere M and its imprisoned m. 
Let V be the velocity of tha/t centre of inertia, II the relative velocity of M and ra. 
If V and R be given in magnitude, R given in direction, V may have any direction, 
and in Maxwell's distribution, for given direction of R, all directions of V are 
equally probable. Conversely, if, whatever be the values of V and R, for given 
direction of R all directions of V are equally probable, Maxwell's law prevails. 
Now consider a very great number of doublets, all having their relative velocity and 
the velocity of centre of inertia within limits R, R + <^R? and V, -V -f dV- Consider 
them before and after collisions between M and m. Nothing is changed by collision 
except the direction of R, and that change of direction is independent of the direction 
of V. Therefore after collision for given direction of R all directions of V are equally 
probable, and therefore Maxwell's distribution prevails after as well as before col- 
lision, and is therefore not affected by collisions. 

To proceed to more general cases. 

2. The characteristic property of collisions of conventional elastic bodies is that 
with continuous variation of the coordinates, and without variation of the kinetic 
energy, the velocities at a certain instant change discontinuously. The general 
treatment adapted to systems of this kind is as follows : — Let there be a system 

defined by n coordinates p 1 . . . p n , the corresponding velocities being p l . . . p n) and 

the generalised components of momentum q l . . • q n . At a certain instant a collision* 

ix.) discontinuity of p x * . , p n occurs. After the collision, let the velocities and 
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• »» * 

momenta be denoted by p\ . . . p' ni q\ . , . q n . Let IT be the potential. In the con- 
figuration Pi • • • p n there are n — 1 independent linear functions of the w forces. 



* • « 



^] 



such that 



&jPri 



, each equal to zero. We might then find n new coordinates c x . . . c„, 






. . . 



^u 



dfe 



»-i 



of momentum corresponding to c 



0, and therefore if S 1? S 2 . . . S* be the components 



* « • "L/ 



are finite or 



zero. In the limit when the remaining force becomes infinite and acts for an infinitely 



short time t 9 



pdS r 



Q (it 

original coordinates. 



dt=0 or S'i — St = . 



1/ 



«-i 



S^^jl = 0, and restoring the 



a i (pi - i>'i) + &i (k -i>\) + • ■ • + h G>f 

» » « * * * 

a„. (p 1 - p\) + . ." . 
in which the coefficients are functions of the coordinates. 



p'.) = <r 

= 



)- . . . (A), 







J 



And therefore n — 1 linear functions of the p's are unaltered by the collision, namely :■ 



<hPi + h \P% + • • • + hPn = %P'l + &]^' a + • • • + *iP» = -! 

«a2>i + &2P2 + • • • + * 2 ^» = a zP'i + ffsPa + • • • + hp'n - S a 

a # •»##*$» « * © # a 9 * e s a * # 

» • # e * * 

a «-iPi + &*-i JPa + * • • + *»— iJP* = a*-ijp'i + ^a-iP's + • • • + kt-iP'n — S 
Again ? since the kinetic energy is unchanged, 



suppose 



**\ 



«-i 



>>. (B). 



Xpq = 2p g , 



and by the properties of generalised coordinates 



Therefore 



* • 

%p'q = %pq 



S (? + ?') (i> - jV) = 



8 t B . ., * . * • 



(C). 



The last equation forms, with the n — 1 equations (A), a system of n equations, 

all linear as regards _p 1 — p' l9 &c. Since p x — p\> &c, are not all zero, we must 
equate the determinant of the system (A) and (C) to zero. That gives us a linear 
equation in q x + q\ } q% + <£& &c- We can now substitute for the q's their values in 

terms of p x • » . jp„, and so obtain a linear equation. 
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a (Pi + Pi) + P(P*+ P'*)> &c - = °> 
or 

a£ x + £p 2 + ...== K == - (ap^ + /3p' 2 . . .), 

where a, /?, &c. are functions of the coordinates and constants. 

That is, we have in all n — 1 linear functions of the velocities, namely, Sj_ . . • S*_!, 
which remain unaltered by the collision, and one other linear function, R, which 
remains unaltered in value, but changes sign. That must be the case on every 
collision of elastic bodies. 

3. The kinetic energy may be expressed in terms of the n variables S x . . . S w _ 1? 

and R, in lieu of the n velocities p x . . . p n , and since it is not altered by the collision, 
which changes the sign of R, leaving S x . . . $ n -i unaltered, it must be of the form 

2E=/(S 1 ...S„_ 1 ) + XK*, 

in which /(Sj . . . 8 nmmi ) is a quadratic function of S 1 . . . S»_i with coefficients functions 
of the coordinates and constants of the system, and X is a function of the coordinates 
and constants. 

4. The system, after collision, has velocities p\, &c, which we will call the 
second state. We may conceive a system with the same coordinates having 

velocities — p\, —p\> &c., and call this the second state with reversed velocities. In 
this state, Sj . . . & n .~i will have opposite signs to those they have in the first state, 
and R has the same sign as in the first state. The system retraces its course, and a 

collision occurs changing — p\ into — jp 2 , &c., leaving S x . . . 8 n ^i unaltered, and 
changing R into — R. 

5. To define a collision, we may suppose that a certain function, \jj, of the coordinates 
and constants is prevented by the physical conditions of the system from becoming 
positive. When xjj becomes zero, dxp/dt from being positive becomes discontinuously 
negative, and a collision is said to take place. We may take xjj for one of our 

generalised coordinates in lieu of p tt9 and \jt, or d\p/dt, for the corresponding component of 

velocity. The kinetic energy is a function ofp 1 . . . p n -u ^ 3 and we may express it as 

a function of S x . . . & n _ i9 and \jj, where S x . . . S„_i are the constants found above. 

Since the kinetic energy is not altered by the discontinuous change in xjj, whatever the 

values of S x • . . S«_!, it must be of the form f(8 1 . . . & n _i) + ^X^. That is \j* is 
reversed in sign, but unaltered in magnitude by the collision, and is, therefore, equal 
or proportional to R found above. 

6. What we have proved for a system is of course true if for system we write pair 
of systems. For instance, let there be two sets of systems : (1) systems M defined by 

« • 

coordinates pj . . . p n and velocities p 1 . . . p r , and (2) systems m defined by coordi 

• * . 

nates and velocities p r + 1 . . . p H and p r + l . . , p n . If xb is a function of p l . . . p n such 
mdcccxcii. — a. 3 a 
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that when xjj = a collision — -i.e., a discontinuity of the velocities — occurs, we may- 
treat the pair of systems M, m in all respects as a single system within the preceding 
investigation. 

7. All those systems, or pairs of systems, for which at any instant \f* lies between 
zero and — (dijj/dt) $t, d\jt/dt being positive, will undergo collision within the time St 
after that instant. We may, therefore, take dxjj/dt or R as measuring the frequency 

of collisions for given values of p l9 &c. 

8. From the linear equations (B) above given we can find any of the velocities, for 

instance, p }} as a linear function of S x . . . S w _i R ? and p\ will be the same function 
with — R written for R. Hence p 2 — p\* = 4R2ftS where 

S/xS == /x 1 S ] + /x 2 S s + &c, 

and the pa are functions of the coordinates and constants. 

Also (p* - p\*) R = 4R 3 %S. 

Now, without altering E or R, or the coordinates, let us make S 2 . . , S«_i pass 
through the whole range of values consistent with 

2E = X 3 R 3 +/(S 18 . . S^O (E). 

Also let <j) (Si . . . S w _i) dS± . . . cZS M _i be the number of systems for which these 
variables lie between the limits 

S]_ and Sj + ^S, 



b^„i and b^_i "|- do 



H—U 



E and R and the coordinates being constant. 

X ULt?U 

• * • {Pi ~~~ ¥1) R . ^ (S x . . . S»-\) dSj . . . dS, 

. . . <j& (S x . . . Sfl_i) S/xS cZS x . . . dS^_i, 



«-i 



4R* 



the integrations being over all values consistent with (E). 

Now, in the Maxwell-Boltzmann distribution <j> (S 2 • . . S^) is a function of 
the kinetic energy only, and is therefore constant throughout this integration. There- 
fore the right-hand member of the last equation is zero, because for any given set of 
values of S } . . . S„_i satisfying (E), there is a corresponding set with the opposite 
signs also satisfying it. Therefore 

. . . R(£ (S-l . . . S*_i) (p x 2 — Pi) dS l . . . dS«_i = 0, 

and, therefore, the average value ofp ± 2 — p\ 2 for all collisions given E and R is zero. 

* • 

The same is true for p£ — p'£ 9 &c. 
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And, therefore, since E and R are arbitrary, p^ 9 &c., are not altered by collision 
at all, that is, the Maxwell-Boltzmann distribution,, given existing, is not altered 
by collisions. 

The above proof also shows that p 3 — p\ 2 is zero without the factor R, that is, the 
average value for all systems is zero, as well as for all collisions ; and in proving that 

p% _ pt^ ^ ( j oeg no |. ma tter whether we introduce the factor R or not. 
9. We will give certain examples of the functions S x . . . S» m1 R. 

I. Elastic spheres of masses m and M respectively. Here a colliding pair, which 
corresponds to a system in the general treatment, has six degrees of freedom. There 
should, therefore, be five linear functions of the velocity unaltered by collision. They 
are x, y 9 X, Y, the tangential components of velocity at the instant of collision, and 

wiu + MU = (M + m) V = mu + MU', 

where u 9 U are the normal components. 
We have from the last equation 

m (u — u') + M (U - U / ) — 0, 
and by the equation of energy 

m (^a - u<%) + M (U 2 - U /3 ) = 0, 

whence 

n + u ' = U +U', 
or 

u _ U = - (u - U ; ) = R. 

II. The system consists of a sphere of mass m colliding with a spheroid of mass M. 
It is assumed that the spheroid will acquire no rotation about its axis of figure, but 
may have rotation about any other principal axis. It has then five degrees of 
freedom, and the system of sphere and spheroid has eight. 

We require, then, seven linear functions of the velocities to be invariable. 




Let be the centre of the spheroid, OZ in the plane of the paper, its axis of 
figure, P the point of collision, PN normal at P, ON perpendicular to PN, and 
ON = c, A the moment of inertia of the spheroid round an axis through O perpen- 
dicular to OZ. Then our seven constants are 

3 G 2 
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(1), (2), (3), (4). Four tangential components of velocity. 

(5). The angular velocity 6 of the spheroid round an axis perpendicular to OZ in 
the plane of the paper. 

And if u, IT be normal components of velocity, and oj the angular velocity round an 
axis through perpendicular to the plane of the paper, the following two, viz. : — 

By conservation of momentum, 

(6.) mu + MIT = (M + m) V = mn + MIT'. 

By conservation of moment of momentum round the axis through perpendicular 
to the plane of the paper, 

(7.) - McU + Aa> = S = - McU' + Aa/. 

We then form the equations 

m(ur-u') + 'M{\J-'U') =0, .- 

- Mc(U - 17) + A(co - a,') = 0, 

and by conservation of energy, 

m (u + ^0 (w ■—■ ?0 + M (U + U') (U — XT) + A (co + g>') (w — a/) = ; 
and equating the determinant to zero we obtain, neglecting common factors, 



u _ y —CO) 



^ — U — Co) 



and, using V, S, and p for S x . . . S w-1 and R in the general equations, we see that if 
the Maxwell-Boltzmann distribution of energy exist, it is not disturbed by colli- 
sions between spheres and spheroids. 

III. Professor Btjrnside's problem (see his paper, Roy. Soc, Edinburgh, July 18, 
1887). He supposes a number of similar and equal spheres, each of unit mass, but 
each sphere, instead of being homogeneous, has its centre of inertia at a distance c 
from its centre, c being supposed very small compared with the radius. The principal 
moments of inertia are for each sphere A, B, O, and the direction cosines of c referred 
to the principal axes through the centre of inertia are for each sphere the same, viz., 
a, yS, y. The direction cosines of the line of centres at impact referred to the principal 
axes are for one sphere L, M, N, and for the other Z, m, n. 

Further the normal velocities are IT, u 9 and the angular velocities round the 
principal axes are Q, l9 fl 3 , Xl 3 , co v w 2 , o> 3 for the two spheres respectively. 

Finally we write, 

N/3 — My = P n/3 — my = p. 

Ly — Na = Q ly —< na = q. 

Ma — L/3 = R ma — Z/3 = r 
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The system of two spheres has twelve degrees of freedom. We require, there- 
fore, eleven linear functions of the velocities to he invariable. They are as follows, 
viz., four components of velocity in the common tangent plane, x, y, X, y, and seven 
others, viz., 

. u + U = V =u f + U', 



cpu 


+ Aco x 


= S 1 


= opu 


+ Aco^, 


cqu 


+ b<j02 


= % 


= cgw' 


+ Bc/ S , 


cru 


+ Ca> 3 


= h 


= m/ 


+ Co/ 3 , 


cPU 


""■" xLUti ■ 


= ®1 


^cPU' 


~~" .£1l<I2 -j, 


cQU 


- BI2 3 : 


= S 3 : 


= cQU' 


— Ba' 2 , 


cEU 


- CO3 : 


= S 3 : 


= cRIT 


~~" OiJ/ g. 



As before, we form the equations, 

w _ tt ' + U - U' =0, 

cp (^ — t^) + A (a> x — o) x ) =0, 

c<£ (it — w') + B (o> 2 — o/ 2 ) = 0, 

cr (a — it') + C (co 3 — o/ 3 ) = 0, 

cP (U - U') - A (% - X^) = 0, 

cQ(U-U') - B (fl a - X2' 2 ) =0, 

cR (U - U') - C (0 3 - 0' 3 ) == 0, 

and by the conservation of energy, 

(u + w ') (^ _ tt ') + (U + U 7 ) (U — XT) + A (o> 1 + a/i) («! — (o\) + &c. = 
and equating the determinant to zero, 

u — U — c (pc^ + go> 2 + ro) 3 + Pfli + QO s + R0 3 ) = />, 
w ' _ XJ 7 — c (po)\ + qct) f Q + rG>' 3 + PO' x + Qf2' 2 + RO' 3 ) = — p 

We can now substitute V, s x> s 2 , s 3i S x , S 3 , S 3 , /> for S x . . . S W-1 R in the general 
equation, and we obtain, as before, the result that the Maxwell-Boltzmann distri- 
bution, given existing, is not affected by collisions. 

10. Professor Burnside obtains the same eight equations as above given, and I 
acknowledge my obligation to him, but he originally deduced the result that the 
energy of rotation is twice the energy of translation, instead of equal to it, as, 
according to the theory, it should be. He has since seen reason to change his views 
with regard to this problem. 
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The equations given by Professor Burnside can easily be modified so as to meet the 
case of elastic bodies of any shape. 

11. It would not be difficult to extend the method of (8) and show that the 
Maxwell-Boltzmann distribution is a necessary, as well as a sufficient condition for 
stationary motion. But that is more completely done by following or extending 
Boltzmann's method. 

Let there be a set of systems which we will call systems M, whose co-ordinates and 

velocities are jp-, . . . p r and p l . . . p r . 

Let the number of such systems for which at any instant the co-ordinates lie 
between the limits 



p l and p l + dpi 



> 



p r and p r + dp r 



$ • 



• # * 



CESJ? 



and the velocities between the limits 



» • s 

p l and pi + dp 



• • » • B fi 



p r and p r + dp r 



be 



or, shortly, 



F (p Y . . . p r p } . . . p r ) dp l . . . dp r dp l . . . dp n 



¥dp l . . . dp r dp l . . . dp,„ 



(A), 






Let there be another set of systems, which we will call systems m, whose co- 

_ • 

ordinates are p r+1 . . . p n and velocities p r + 1 . . . _p«. 

And let the number of systems m for which at any instant the coordinates lie 
between the limits 

p r + i andp f + 1 + dp r + 1 



e * i * 9 



p u and p n + d!p w 



* e 



r « « * • . o . »\ _L_J /. 



and the velocities between the limits 



~\ 



p r+1 andjp r + 1 + dp r + l 

* B * • • * 9 9 f 1 ^ 

9 4* » 

p H and p n + dp,, 



a o » a ft 



(B'), 



be 



/ \pr-\-\ * • • Pn Jpr + 1 * • » jPw/ ^JPr + l • • • ^p« tyP r + 1 * • • ^P?» 
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or, shortly. 

The number of pairs of systems, each consisting of one system from each set, whose 
coordinates and velocities at any instant lie within the above limits, is 

dp x . . . dp r dp r+ i . . . dp fl Ff. 

Now let iff be a function of the coordinates which cannot become positive, and 
sach that when xjj = the velocities change discontinuously, and a collision occurs. 
We may use xjj for one of our coordinates, expressing p n in terms of p x . . . p n _ x 

and \jt, and in like manner p n in terms of p x . . . p^_! and i//. All those pairs for which 

at the given instant xjj lies between zero and — \p8t, \p being positive, will undergo 
collision within the time St after that instant. And so the number of such collisions 

in unit time is, writing R for \ft, 

• * _ _ 

dp x . . . dp n _ l dp x . . . dp n _ :i dR . F/R. 

Here F is a function of p x . . . p r only, but/, by virtue of the elimination of p n and 

• » • 

p m is now a function p x . . . p nm .\ and^ . . . p n ~\ R- 

« 

Each of these collisions changes p x intoi 9 ^, &c, that is, changes the system from 
the first into the second state. The number in unit time of collisions, which with the 
same coordinates change the system from the second state with reversed velocities 
into the first state with reversed velocities, is 

dp x . . . dp n _ x dp\ . . . dp' n _i cfftF'/'R, 

in which F', f are the same functions of p\, &c, as F, /.are ofp 1? &c. 
By a known theorem 



• • • 

(X u j . • . (X u yi j —— Cva/i ... 



n — 1 ? 



and in Maxwell's distribution F/'= F/'. And so the number of reverse collisions 
is in that distribution equal to the number of direct collisions. And this insures the 
permanence of the distribution. It is assumed that there are always as many systems 
with any given set of velocities as with those velocities reversed. And so we speak 
of the second state with reversed velocities as equivalent to the second state. 

12. If F/' =£ I/, the number of reverse collisions is not equal to the number of 
direct collisions. And, therefore, more (or fewer) pairs of systems pass out of the 
first state into the second than vice versa. In that case the number of systems M, 
whose coordinates and velocities are within the limits A, A 7 , is increased by collisions 
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with systems m 9 whose coordinates and velocities lie within the limits B, B', by the 
quantity 

dp 1 . . . dp n -i dp 1 . . . dpn-\ dR (Ff — Ff) R per unit of time, 

and is increased per unit of time by collisions with systems m for all values of 
Pr + 1 • • • Pn by the quantity 

dp 1 . . . dp r dp i . . . dp r . . . (Ff — Ff) R c?p,. + 1 . . . dp n _ l dR, 



all values of p r + u &c, being included in the integration. 

We will now assume (see p. 418, post) that the velocities of M and m are not on 
average altered except by collision between M and m. And so the above-mentioned 
increments are the only increments by which the class of M systems within the limits 
A, A' is affected. 

In that case 



dF 

dl 



dp x . . . dp r = dp x . . . dp r l . . . (Ff — Ff) R <%>?<+ 1 . • • dp nm _ t dR } 



and, therefore, 



• • • 



dF 

clt 



log F d^9 T . . . dp r (over all values of p Y . . . p r ) 



. . . (Ff — Ff) R log F dp x . . . dp nmm i dR. 



*/ w 



By symmetry, as the right-hand member includes all possible collision between 
and m 9 

df 

. . . - d jlogfdp r+l . . . dp* 



• .• • o^ -" *oo -^ i°g/ c ¥i • • • ^»-i ^> 



and, therefore, 



»' w 



• » • 






log F dpj . , . dp r + . . . j. l°gfdp r +L . . . dp n 



. . . (Ff — F/) R log (Ff) dp x . . . dp w-J dR. 



By symmetry, as we may interchange the accents, 



dF 
dt 



» 

log F dp x . . . dp f + 



8S» 



^ 

* 



log/rfp f + , . . . dp H 



IT . . . (F/- F/') E log (F/') dpi . . . dp„Li dE, 
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and. therefore, 



i 

2 



w 



j j j . . . (F/ 7 — F/) R log ~/~ dp l . , . cVp nmml dR: 



Now let 



H = . . . F (log F — 1) dp l . . . dp r +|| • • • /(log/ — 1) dpr+i • • • dp«> 



and, therefore, 



dt 



[[[ dF [ 



r f 



« » • 






logfdp r + 1 . . . dp n 



l 

2 



w 



. . . (F'/' — F/) R log ^7, #i . . . dp«-i ^R, 



as we have seen. 

Now, this expression is necessarily negative, unless F/ 7 = F/ whenever the pair of 

systems having coordinates and velocities p x . . . p n can pass by collision, and, there- 

fore, with unchanged kinetic energy, into the state in which they are p x . . . p' n , that 
is, unless the Maxwell-Boltzmann distribution exist, and is then zero. H therefore 
tends to a minimum which it reaches when F/ = F/ 7 . 
We will now make 

H = H-j_ + K 



where Hj^ is the minimum value assumed by H when F/' = F/! Then dH/dt = 0, 
and dK/dt = dHjdt, and is always negative. We may define the function K to be the 
disturbance, and (1/K) (dK/dt) to be the rate of subsidence of the disturbance by 
collision. In certain cases, we can calculate the rate. 

13. We have assumed that / varies only as the result of collisions. That is, if 
df/dt denote the time variation of / due to causes other than collisions, and dHjdt be 
formed from df/dt as dHjdt from df/dt> then 3H/9f = 0, on average. It is worth 
while to consider on what condition this may be safely assumed. 

Let 



3H 

"dt 



tif 



As we are dealing with rigid elastic bodies under the action of no forces, we may 
treat / as a function of three translation velocities, and three angular velocities, 
iv l9 iv 2 , Wfy about three principal axes of the solids. Let A, B, C, be the principal 
moments of inertia. Evidently, there being no forces, the translation velocities cannot 
vary except as the result of collisions. But for each solid, iv l9 iv 2 , w% may vary, the 
law of their variations being Eulek/s equations. We may, therefore, in calculating 
dH/dt treat /as a function of iv i9 iv 2 , tv 3 only. Then 

MDCCCXCII. — A. 



3 H 
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3/ ^^ df dw l df dw 2 df dio« 
dt dw Y dt d%o% dt dw s dt 

B-C df t C-A df , A - B df 

1 A o 



and 



9H _ B - C 

¥ — ~ A ~ 



l~— log/, W 3 W 3 . tfw x C?W 3 ^3> 



with two other corresponding terms, the limits being in each i oo . Now, with these 
limits, 



■j— log f. w % w B . dw x dw 2 cliOo 



j <% dw s iv 2 w 3 { f log /,, = co - / log / Wl =: « „ } . 



Now, we may assume /*= and /' log /* = 0, when any one of the three variables is 
infinite, whether positive or negative. And this assumption is sufficient to justify 
the statement dH/dt = 0. 

14. It is possible to calculate, in a simple case, the rate at which a disturbance subsides 
by collisions. For example, two sets of elastic spheres, N of mass M, and n of mass 
m, in unit of volume. In the normal state, the number in unit of volume, whose 
velocities are represented by lines drawn from an origin to points within an element 
of volume U 3 sin d6 d<j> dU is for the M spheres 

AM\f 

NF e~ mija Wmn0d0d<f>dTJ, 
where U, 0, <j> are usual spherical coordinates ; or, let us say, 



F(U) = N^Ye-* MD * 



Similarly, for the m spheres, 






expresses the law of distribution of velocities in the undisturbed state. We will 
write F and f for these expressions. 

We will now suppose there is a small disturbance consisting in h having different 
values for the two sets. Let h be written h (1 + D) for the M spheres, and h (1 + d) 
for the m spheres. We shall neglect third and higher powers of D, d. Then F 
becomes in the disturbed state 



(fp + D > ! 
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that is. 



N( — J (1 + J)fe- m{] "{l 



hDMXJ 2 + terms in D 2 }. 



(It will appear that terms in D 2 ? &c, are not required.) 
Similarly /becomes 



n 



Ivm 



IT 



(1 + df e - *»«*(! — ctawm* + eP &c). 



We will further suppose that the disturbance is introduced without changing the 
total energy. That gives the relation 



N. 



2A (1 + D) ^ 2h (1 + D) 97 > V T /? 



2A 



or 



N 



•■•Im 



'/t 



1+D 1 + d 



N + w- 



15. The disturbance will subside by collisions between M and m. And we will 
treat of the case in which it subsides in such manner that the above values of F and/ 
apply at every instant with the values that D and d at that instant have. Such a 
mode of subsidence is possible, at all events if our equations lead (as they do) to 
a relation of the form (1/K) (clK/dt) = constant. 

Let us then form the function 






J J 







. F (log F - 1) IP sin a da dp dU, 



where U, a, and ft are usual spherical coordinates, 



'oo ,"7r C2tt 



+ I /(log/ — 1 ) v? sin a da dfi du, 

Jo J o J 



where F and /have the values above given. 
That is 

H = f log (1 + D)f f fF.lP sin a da dp dU 

j j j 



or 



because 



+ f log (1 + d)\\ f.ur sin a da d/3 du 
-f- terms independent of D and d, 

H = N|log(l + D) + rcflog(l +d), 
-j- terms independent of D and d, 



FW sin a da d/3 dU =N 



«/ W » 



/^ 2 sin a c/a a/3 d/w 
3 II 2 



n. 
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The terms independent of D and d are equivalent to H 1? the minimum value of H 
when there is no disturbance. 
And so 

K = H-H 1 = fNlog(l + D) + fnlog(l + d). 



That is 



K = f JcZ-f) + -§N(l) 






Now since 



_j_ — - r= N + ^, 



1 4- D ' 1 •+ d 

ncl 



D 



N" + (NT + w) d! 



and 



K — -§-« /V/ — -A — £ N - 7 - — - — £ — 

2 \ 2/ a N + (N + w) d * ' {N + (N + rc) cZ} 3 

= £ £ (N + n) d\ 

In order to find dK/dt we will transform our coordinates thus : Let V denote the 
velocity of the centre of inertia of a pair of spheres M and m, p their relative velocity, 
9 the angle between V and p. Then 

1 \M + 7/i W M + m r 



and 



and so 



MIP + ina* = (M + m) V 3 + ( - — 



.2 



F/= N« f— Y f— y (1 + D)* (1 + d) : U- h i^ + m>T2 + w^} (1 - DAMU 8 - dhmu*) 



F'f = N»i (^Y (^Y (l + D)* (1 + d)» e" A (< M + w)V2+ ari7/} (1 - DM1U' 8 - &»' 2 ) 
where U', it' are the values of U, u after collision. 
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Therefore 

Vf-Tf = Nn(^^)*(l + D)»(l +^ f 6-^ + *> Vi+ 5^^ 



{DAM(U'* - IT 3 ) + dAm («* -«»)} 



Now, if 0' be what (9 becomes after collision, 



U' 3 — IP = 2Vp . — — — (cos ff — cos 0) 

M + m x } 



it* ~~ y? ; = ~~ 2Yp r-p^ — (cos & — cos $) 



Also by the relation 



1~v / vtv 



N + (N + ») ^ ' 



and making these substitutions 

F/' - F/ = Nn(^'(^-)*(l + D)«(l + d)«€-*{<"*-^-M^} 

' Mm Tr N + n , \ 

A \iT7~7 2 • V P • tSTTxTT^Tx'V ( cos " — cos 0) a . 

Also 

1 Ff Mm N + 7* , A nn 7 

lOg =77? = A — 2 Vp- T7T ~ ri — (COS — COS ) cZ. 

In forming (F/' — F/) log (Ff/F'f) we see that the last factor is squared, and so 
the product contains the factor <p. We may, therefore, now write 1 for (1 + D) 1 and 
(1 + df, and also write (N + 7i?)/N for (N-.+ n)/{'N + (N + n)d}, otherwise we 
should have terms in #. 

Therefore 

(F'f - F/)log^> = - Nn(^J(^ J € -'{« + ^ + ^} 



A * * (i^T 4V V "4^ (cos 0' - cos *)». 



M + m/ r N 



Again, d'K/dt contains the factor \fj denoting the frequency of collision. Now the 
number of collisions given V and p in unit of volume and time is proportional to 
7rs 2 p, where s is the sum of the radii of M and m. Also all directions of p before 
collision are equally probable, and given the direction before collision, all directions 
after collision are equally probable. Therefore, given V and />, the number in unit of 
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volume and time of colliding pairs for which the angle between V and p before 
collision is between 9 and 9 + d9 is \irs 2 p sin 9 d0, and the number for which after 
collision it lies between f and ff + d9" is ^rrs 2 p sin & d9\ Hence we have to 
multiply (cos & ~ cos 9) % by ^rrs 2 p sin 9 d9 sin & dd\ and integrate in each case from 
ir to 0. The result is f ns 2 p. 
And so we get 



W \M + m) \ IT J \ 7T / 



cTK. ^ 2 
~dt ~ ~ ¥ 

'27T r°° c T „ M*« 

^ rfv sm <x da dp 
' ->o 





/*QO 


fTT 


X 






fc 


" 


o J 



V s d V sin a da d{3 dp 4tt P 3 e~ '' l ( " + "' Vi + a " + ""* } 4 V 9 p 8 . ™ 2 



3 



where V, a, /3 are usual spherical coordinates 

Also, as we have seen 

and therefore 

K d*' — 3\/7r ^ + U ' (MTm) 1 a/7* 

— — C suppose. 

And if K , D , d be initial values, 

K=Ko€- c ', D = D €-* C S d = d €~* ct , 

the rate of subsidence is directly proportional to the density and to the square root of 
the absolute temperature.* 



* Since writing the above I find that this result has already been obtained for the case of elastic 
spheres by Professor Tatt, by an independent method. 



